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, Abstract. We consider integrable systems that are connected with orthogonal 

separation of variables in complex Riemannian spaces of constant curvature. An iso- 
morphism with the hyperbolic Gaudin magnet, previously pointed out by one of us, 
' extends to coordinates of this type. The complete classification of these separable 

^ , coordinate systems is provided by means of the corresponding L-matrices for the 

Gaudin magnet. The limiting procedures (or e calculus) which relate various degen- 
erate orthogonal coordinate systems play a crucial result in the classification of all 
such systems. 



PACS: 02.20.+b, 02.40.+m, 02.90.+p, 03.65. Fd 

1. Classical Integrable Systems on Complex Constant Curvature Spaces 
and the Complex Gaudin Magnet. Separation of variables in the Hamilton 
Jacobi equation 



■A ^ dW 
(1.1) H{pi,...,pn;xi,...,Xn) = g"^Pc^p/3 = E, Pa = ^ — . a = l,...,n 
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amounts to looking for a solution of the form 

n 

(1.2) W =Y,Wo.{Xa;hi,...,hn), hn = E 

a=l 

The solution is said to be a complete integral if (iet{d'^W / dxidhj)nxn 7^ 0. The 
solution then describes that of free motion on the corresponding Riemannian space 
with contravariant metric gr"'^. Indeed, if we require hj = — hntSnj, j — 1, ...,n 

for parameters bj, we find that the functions x(b, h), p(b, h) satisfy Hamilton's 
equations 

(1-3) Xcc = Pa = --^ 

OPa OXa 

In this article we allow the Riemannian space to be complex and we consider variable 
separation of (1.1) on the following two classes of spaces. 

(1) The n dimensional complex sphere SnC- This is commonly realised by the 
set of complex vectors x = {xi, ...,Xn+i) which satisfy X^^ii = 1 and 
have infinitesmal distance dx ■ dx = y^"^i dx"^ 

(2) The n dimensional complex Euclidean space Enc- This is the set of complex 
vectors x = (a;i, Xn) with infinitesmal distance cZx • dx = dx"^. 

A fundamental problem from the point of view of separation of variables on these 
manifolds is to find all "inequivalent" coordinate systems. As yet, this is an un- 
solved problem, principally because many such coordinate systems are intrinsically 
nonorthogonal. For orthogonal coordinate systems the problem is completely solved 
and in this case the constants hi occuring in the complete integral can be chosen 
to be the values of an involutive set of constants of motion 

n 

(1.4) Aj^ ^ aalPaPp, j = 'i-,-,n, An = H, 

{Aj,Ak} = 0, 

where 



Si7( \r( M X-fdF dG dF dG 
{F{x^,P,)Mxa,P,)} = - a^a^)' «'/^ = 1' 

is the Poisson bracket. These constants of the motion are such that 

(1) each of the tensors a^^ is a Killing tensor and satisfies Killing's equations 

V(a4'!)) = 0, [1], and 

(2) the Aj can be represented as a sum of quadratic elements of the enveloping 
algebra of the Lie algebra of symmetries of each of these two considered 

spaces. 

The Lie algebras of these spaces have respectively bases of the form 

(1) SO{n + 1) : Mo,p = Xapp - xppc, a, /3 = 1, n + 1. Here, 

(2) E{n) : M^p, P^=p^, a, 7 = 1, • • • , n, a 7^ Here, 

{Map, Py} = Sp^Pa - SayPp, {P«, P/l} = 0. 
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In this section the separable coordinate systems classified in [1,5,10] arc given 
an algebraic interpretation. This is done using the complex analogue of the isomor- 
phism between all integrable systems connected with all possible separable systems 
and the m site 5*0(2, 1) Gaudin magnet, [2,4]. The m-site complex Gaudin magnet 
can be realised as follows. Consider the direct sum of Lie algebras, each of rank 1, 

(1.5) ^=e-^iSo«(3,C). 

The generators Sq, e C^, a = 1, • • • , m of ^ satisfy the Poisson bracket relations 



(1-6) {si,4} = -Sapeijkst 

The following metric will be used subsequently when norms and scalar products 
are calculated: 

si = {Sa, Sa) = (sl^f + {slf + {slf, (s^, S^) = S^sJ + 4^ + S^S^ 

If for each a, = then the variables lie on the direct product of n complex 
spheres in C^. The complex Gaudin magnet is the integrable Hamiltonian system 
described by the n integrals of motion Hq, which are in involution with respect to 
the Poisson bracket, 

(1.7) H^ = 2J2^-^^, {H^,H,} = 0. 

(We will give a simple proof of this involution property later.) Here the 6^ are 
taken to be pairwise distinct. This integrable Hamiltonian system is called an m- 
site S0{3,C)—XXX Gaudin magnet. The are all quadratic functions in the 
generators of the A algebra and the following identities are satisfied: 

mm m 

(1.8) J]iy« = o, ^e«iya = J'- 

a=l a=l a=l 

where we have introduced the variables 

m 

(1.9) J=E^«' J' = (J'J)' 

a=l 

the total sum of the momenta Sq. Indeed 

(1.10) {J^J^} = -eyfeJ^ {j\Ha} = 0, i,j,A; = 1,2,3, a = l,...,m. 

The complete set of involutive integrals of motion is provided by Ha, and, for 
example, J^. The integrals are generated by the 2 x 2 L matrix [3,4,6,7,8], 
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where 

(1.12) detLM^-(^Vi.C) = -i:^-E^ 

a=l a=l 

Furthermore, L{u) satisfies the hnear r matrix algebra [2,21] 



(1.13) {L{u),L{v)} = [P,L{u) + L{v)l i = ^ 

u — V 

where 



(1.14) P = 



/I 
















1 








1 








Vo 








1/ 



1 2 

L(u) = L{u) ® /, L(w) = / (g) L{v). 



The algebra specified by (1.13) contains all the necessary Hamiltonian structure of 
the problem in question. Note that (1.13) is equivalent to the easily proved relations 

(1.15) 

{A{u),A{v)} = {B{u),B{v)} = {Ciu), C{v)} = 0, 

{A{u),B{v)} = ^—{B{v) - B{u)), {A{u), C{v)} = ^—{C{u) - C{v)), 
u — V u — V 

{B{u),C{v)}^-^(A{v)-A{u)). 

Lb U 

From (1.15) and the Leibnitz property of the Poisson bracket it is straightforward 
to deduce that 

{detL(M),detL(v)} = 0. 

In his article [4] Kuznetsov has explicitly given the nature of the isomorphism 
between the XXX Gaudin magnet models and the separation of variables on the 
n dimensional real sphere Sn- The purpose of this article is to extend these ideas 
to complex orthogonal coordinate systems on the complex n sphere Snc and, of 
course, as a consequence complex Euclidean space Enc Following Kuznetsov [4,7] 
in the case of the sphere, we set Cq, = 0, a = 1, • • • , n + 1. The coordinates on the 
resulting cones are parametrised by 

(1.16) 4 = \{pI + ^l)^ 4 = \iPl- ^a)> 4 = l^PocXoc- 

It follows from (1.6) that {xq., xp} = {pa,P[}} = 0, {pa, xp} = dap- Introducing the 
new variables Mq,^ = x^Pp — xppa which are the generators of rotations we have 

These generators satisfy the commutations relations given previously. This equality 
establishes the simple quadratic connection between the generators of A and the 



SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 



5 



subsequent discussion) transform into the following integrals for the free motion on 
the n sphere 

Ha - rip 2 



For = we obtain the Casimir element of the so{n +1) algebra X^q;<^ 
The total momentum J takes the form 

= ^(p-P + x-x), J2 = -(p-p-x-x), j3 = -p-x 

where the scalar product for the vectors x and p in C"'"'"^ is Euclidean. The quantities 
M and J form the direct sum so{n + 1) ® so(3) as a result of the commutation 
relations 

{M,^,j^} = o, {J^J^■} = -eyfcJ^ 

Therefore, in addition to the involutive set of integrals i^Q, we can choose = 
i Yla<i3 ^l*^^ + ^'^^) = X • X = c, which givcs the equation of the n sphere. 

2. Generic Ellipsoidal Coordinates on Snc and E^c Critical to the sep- 
aration of variables on the n sphere S^c is the system of ellipsoidal coordinates 
graphically pictured by the irreducible block 

(2.1) ('S'„c|ei|e2| ■ ■ ■ |e„+i|. 

where in general ^ for a ^ (5. The separation variables are defined as zeros 
of the off diagonal element B{u) of the L matrix, i.e., B{uj) = 0, j = 1, ...,n. It 
follows that 

(2.2) y =0 for u = Uj and xl = c j, ^ ^. 

^^^i-ea n^^c.(e/3 -e^) 

Each vector of momentum Sq, is associated with a cell Ca of the block. Note that 
x^ = 2(s^ + is'^). For each Uj the conjugate variable Vj is defined according to 

n+l 



(2.3) Vj = -iA{uj) = ^ E ;r 



2 „ 1 
oc=l •' 



Prom (1.15) one can show that the Wj, Vi satisfy the canonical relations [4]: 
(2.4) = = 0, {vj,Ui} = bij. 

The change to the new variables Vj,Ui,c and is effectively the procedure of 
variable separation of (1.1) in ellipsoidal coordinates on the n-sphere. Writing the 
L matrix in terms of the new variables we obtain 

fA{u) B{u)\ c n^=^iu-u,) 
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(2.5) A{u)---B{u) ^iJ Ui^j{uj-Ui) 

where A{uj) = ivj, j = 1, ■ ■ ■ , n, and A{u) — >• {l/u)J^ + • • • as u — >• oo. 

To obtain C{u) we first notice that equating residues at ea on the right and left 
hand side of A{u) gives 

This together with the expression (2.2) for in terms of ttj gives C{u) in the new 
variables. Three other useful formulae are 

where the hat in (2.6) means that the product terms with 'j — a and 'j — (3 are 
omitted, 



(2.7) = - E ^ frr'^' \ = E E . 



2 - n-ij(«,-e^) 



These relations together with (2.2) establish the explicit connection between the 
two sets of 2n + 2 variables pa, and w^, Vj, c, J^. The equation for the eigenvalue 
curve r : det(L(it) — = has the form 

-A^ - A{uf - B{u)C{u) = 0. 

If we put u = Uj into this equation then A = ±Vj. Thus variables uj and Vj lie on 
the curve F: 

n 

(2.8) + E _ = - det L{uj) = 0. 

a=l 

Equations (2.8) are the separation equations for each of the n degrees of freedom 
connected with the values of the integrals Hq,- For the sphere Cq, = these have 
the form 
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The Hamilton Jacobi equation (1.1) when parametrised by these variables has the 
form 

which can be solved by the separation of variables ansatz 



(2.10) 



n n „ 

7 = 1 7 = 1 



It will be convenient to employ an alternative form of the L matrix. If we use 
the vector 



n+l 

Sr 



(2.11) L{u) = , L{u) = (s,L(«)), detL(«) = -L{u)-L{u) 

a=l ^ ~ ^« 

where 



-1 or V z or V -1 

we see that L satisfies 



(2.12) {V{u),L\v)} = ^ {L\u) - L\v)) 



At this point we must consider a crucial difference between the real sphere and 
its complex counterpart. In the case of the complex sphere the generic ellipsoidal 
coordinates can admit multiply degenerate forms: the restriction Ca ^ e/3, for ct ^ /3 
can be lifted. The resulting coordinates can be denoted by the block form 

('5nc|e^ \el' I • • • |ej« I, Ai + • • • + = n + 1, 

where the Aq, denote the multiplicities of the e^- To understand how the previous 
analysis applies to these types of coordinates we first illustrate with an example 
corresponding to the coordinates with diagram 

('S'ncleilesI • • • |en+i|. 

In this case we write 

n+l 

(2.13) uu) = + + y 

a=3 

Putting 

(2.14) X2 x[ + 6X2, Xi x[, P2 p'l + ep2, Pi p[, 
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then in the hmit as e — > we find 

Zl 



+ 



Z2 



(u-ei)2 u-ei ^w-ea 



n+l 



where 



(2.15) 



Z2 =< ^(p'lP2 +^'1^2)5 ^(p'lP2 - ^'1^2): +pWi) > 



a = 3, • • • , n + 1. 



The components of zi,Z2 satisfy the £'(3, C) algebra relations 

where the Poisson bracket (expressed in the primed coordinates) is 

""^^ dF dG dF dG 



{F, G} ^^^^dp', dx'^ dx', dp'^ ^ 



j,k=i 



and 



B = (Bjk) 



/O 1 
1 
1 








\0 ••• 1/ 



It can easily be verified from these relations that relations (2.12) are again satisfied. 
Thus, {L ■ L(m), L ■ L(w)} = so the coefficients of the various powers {u — ej)~^ 
in the expansion of L ■ L(u) form an involutive set of integrals of motion. 
We obtain 



L L = 



1 



-Zl • Zl + 



2Z2 ■ Zq, 



Zl 



1 V- 

• Z2 + 7 r^Z2 • Z2 + > ^ 

[u — eiY ^-^ 



2zi ■ Zc 



1 



a>3 
Za,Z/3) 



(w - ei)2(w - Cq,) 



(2.16) 



-1 



+ 



4(w - ei)2 
1 



2(w - ei)(« - Cq) 



n+l ^ 

(p'l^ -p^a;'i)2 + ^ -{j^^x'^ -P'a^f 

^ A{u-eiY{u-ea) 



1 



a,/3=3, 0:5^/3 

In particular 
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To relate this to the projective coordinates on the complex n sphere we recall that 
under the transformation (2.14) the fundamental quadratic forms X = a\x\-\-a2X2-\- 
ES3 P = ^fi + «2Pi + Vl transform to X = 2x'^x'^ + E^ia ^'l. P = 

2piP2 + Z^aisP'a- Therefore if we take the coordinates 

x\ = -^{Xi + iX2), X2 = -^i^^ - ^-^2), x'g^ = Xa, a = 3, • • • ,n+ 1 

P'i = ^iPi + iP2), p'2 = ^{Pi-iP2), Pa = Pa, « = 3, • • • , n + 1 
we then can write 



pW^ - pWi = V2Mi2, P^Y^P^ 



n+1 n+1 

2 

i 

3=1 3=1 



where M,/, = XjPf^ — X^Pj. 

The integrals of motion i^Q,, Z in this case have, using partial fractions, the form 

11 n+1 



n+1 ^ ^ 1 



a=3 

^ 4(ea - e/3) - ei) f::^ 4(ei - Ca) 2 



2(ei - Cq) 

Z Y 



— + 7 + 



^^gW-Ca (tt-ei)2 tt-ei 

where >^ = Ea=3-f^a- 

The analysis presented so far could have been deduced from Kuznetsov [4] where 
the double root is essentially contained in the s systems of type C on the real 
hyperboloid. Furthermore, the threefold root is contained in Kuznetsov's type D 
systems. The question we now answer is how to use these techniques on the case of 
ellipsoidal coordinates corresponding to multiply degenerate roots. For this we use 
the limiting procedures developed by Kalnins, Miller and Reid [5]. We recall that 
the process of using these limiting procedures amounts to altering the elementary 
divisors of the two quadratic forms 

n+1 2 n+1 
^1 
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Theorem 1. Let ui he the generic ellipsoidal coordinates coordinates on the n 
sphere viz. 

(2.18) a^xl= ^, « = !,•••, n+1 

with corresponding infinitesmal distance 

7 2 \ — T-r / \ idtij\^ 

as = — - > lij^iiUi—Un) — XT 

and coordinate curves 

n+l 2 n+1 

(2.19) E^ = 0. E = 

a=l ■' a=l 

Then the degenerate ellipsoidal coordinates having the infinitesimal distance 



(2.20) ds^ = -\ Y.[Tl,^.{u, - u,)] 



2 



can he ohtained from generic ellipsoidal coordinates via the transformations 
e/ ^ ei j = 1, • • • , iVj, J = 1, • • • ,p 



i=2 

Nj 

(2.21) + 



1=2 

where 

'ei-^\_, = W,^^{'e]_, eU), a/ = l/[n,^,('^ej_i -^^ e^.^)], /c = 1, • • • , iV^ 

and A^i + • • • + ATp = n+l. (We require "^Sq = and take the limit as the ^ e\ — > 
for h = 1, ■ • • , Nj — 1.) In particular, 

V Nj 

(2-22) L{u) = yy 7 ' 

where 

i i i 



SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 



11 



(2-24) {(^/),, {zi')m} = -5jJ'{zj^U-Nj)seir 

and we also have 

(2.25) L . -Liu) = \ YP^'A - ^kPj)«Pn - ^nPin)-r— 



[u - ej)p{u - clY 

where 1 < p < Nj,l < q < Nl,1 < j,m < Nj,l < k^n < Nl, subject to 
j +m — Nj + 2—p,k + n = Nl + 2 — q and the summation extends over the indices 
J, L, p, q, j, m, k, n subject to these restrictions. 

From expressions (2.24) we can verify that relations (2.12) are satisfied, so that 

{L-L{u),L-L{v)} = 0. 

The constants of the motion can be read off from the partial fraction decomposition 
of (2.25). This clearly illustrates the compactness of the r matrix formulation for 
the operators describing the integrable systems examined so far. 

In dealing with the case of Euclidean space E^c the most transparent way to 
proceed is as follows. The generic ellipsoidal coordinates in n dimensional complex 
Euclidean space are given by [1,9] 

(2.26) ^, j,a,P = !,■■■ ,n 



with coordinate curves 



(2.27) ^^^ = 1, u = Uj, j = l,---,n. 

Proceeding in analogy with (2.2)-(2.10) we can obtain the r matrix algebra. The 
corresponding L(u) operator is 



(2.28) L{u) = E ;rr 




a=l 

(Indeed the equation L^{u) + iLP'iu) = is just (2.27). Moreover it is obvious, due 
to the fact that expressions (2.11) satisfy (2.12), that expressions (2.28) also satisfy 
(2.12).) The conjugate variables Vj are defined by Vj = —iL^{uj) and they must 
satisfy the canonical relations (2.4). The integrals of motion are determined 
from 



Ha 



(2.29) l.\u) = E — 

where 



1 ^ 



with Ylia Ha = — P^/4. The separation equations are of the form 



•yj + L^(uj) = 0, j = l,---,n. 
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3. Cyclidic Coordinates. Associated with the separation of variables problem 
for the Hamilton Jacobi equation (1.1) with 7^ is the corresponding E = 
problem. In this case the equation is 

-J^ dW 
(3.1) 2^ 9'^'^PaPf3 = 0, P=o — , a=l,---,n 

a,/3=l 

and we consider only complex Euclidean space. While it is true that all the coordi- 
nate systems discussed for E^c with E ^ will provide a separation of variables of 
this equation, there are coordinates that provide an additive separation of variables 
only when E = 0. This is related to the fact that the E — equation admits a 
conformal symmetry algebra [1,11,13]. 

The most convenient way to proceed is to introduce hyperspherical coordinates 

{Xi,-- ■ ,Xn+2}, 



xi = t'^i^ z^j - 1), X2 = it'^i^ + ^^=+2 = /c = 1, • • • , n, 

i=i j=i 

related to the usual Cartesian coordinates {zi, ■ ■ ■ ,Zn} according to 

n+2 

Zk = Xk+2/{-xi-ix2), k=l,---,n ^ = 0. 

We consider the system of §1 in n + 2 dimensions, where J = and the — 0. 
The general (separable) cyclidic coordinates are specified by 

n+2 ^2 n+2 



Furthermore, 



The Hamilton Jacobi equation is given by 

(3 2) = -a^y (dWy u:;±liu,-e,) ^ ^ 

The quadratic forms Cl and ^ have elementary divisors [11 • • • 1], see [5,11]. It 
is known that the geometry of these fourth order coordinate curves is unchanged 
under birational transformations of the form 



CKCfe + P auj + P ^ aX + j3 



t — „ j: n ] ; 1 
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Now we can mimic the exposition given for the Gaudin magnet integrable systems 
using the hyperspherical coordinates and the Poisson bracket 

\ ^/ \-. ^ — ^ / dF dG dF dG , 
(3.3) = T.^^g^^ - 5^^). 

fc=l 

the Xi,pj now being regarded as independent. The analysis then proceeds much as 
in the construction (2.2)-(2.12), but with n replaced by n+1; indeed the coordinates 
are defined as zeros of the off diagonal element 

n+2 2 

B(u) = y 

subject to the restriction $ = 0. The crucial difference is that J = 0. The expres- 
sions for Pa and Map are altered by a factor cr: 



Pa = 2axa^^ 



.^^ ea-Uj ^i^j{Uj - Ui) 

which together with in terms of Uj gives C{u) in the new variables. Another 
useful formula is 

Map = 2axaxp{ea - ep) l^Vj— _ 

In fact the Poisson bracket {,}h can be identified with the Poisson bracket { , 
} for functions defined in the n dimensional space spanned by ^i,--- ^z^- This 
can readily be seen by noting that F{zi,--- , z^Pz^ ■ ■ ■ ,Pzn) — F{—xs/{xi + 

iX2),--- , -Xn+2 /{xi+iX2),-{xi+iX2 )Pxz + (Pxi + iPx2 ),■■■, -{xi+iX2)Pxr,+2 + 

Xn+2{Pxi +'i'Px2)) from which the equality of the Poisson brackets follows identically. 
The infinitesmal distance for general cyclidic coordinates is 



(3.4) ds^ = - 



1 \-2 { TT \,UUi I 



We denote the coordinates defined by this graph as 

[EnC,E = 0|ei|....|en+2|- 
Other coordinates of this type are those corresponding to the graphs 

('S'pc|ei|....|ep+i| e ('S'qc|/l|----|/q+l| 

where p + q = n + 2. These coordinates are given by 
(3.5) 

^f^l(^j - eg) 

n/3^a(e/3 - Cq,) ' 



(^xl = a,p=l,--- ,p+l, k = l,--- ,p 



„ n^=i(t', -/a) 

^^fc = ?F 77 FT' I,--- ,9+1, k = p + l,--- ,p+q + 2 

J-WaU/3 - /a J 
P+1 q+1 

a=-[J2eiX^ + J2fjx^,]-\ 
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We note here that the the fj are pairwise distinct, for if they were not then for 
any of these quantities which occured with multiphcity more than 1 a birational 
transformation could transform it to oo and hence to a graph corresponding to 

EnC- 

It can happen just as in the case of generic eUipsoidal integrable systems on the 
sphere that some of the in (3.5) are equal to some of the fj also. In this case 
the rules for obtaining the corresponding L matrix are summarised in the following 
theorem. 

Theorem 2. Denote the generic ellipsoidal coordinates by the graph 

[Enc\ei\....\en], Bi^Bj 

and generic cyclidic coordinates by the graph 

{CEnc\ei\ |en+2}, ei^Sj. 

Separable coordinates for the Hamilton Jacobi equation (3.1) corresponding to generic 
graphs with multiplicities 

[£'nc|ei'| • • ■ Icp''], mH hnp = n, 

{CEnc\eT I • • • K"), m + • • • + = n + 2, 

[5] can be obtained via the transformations of Theorem 1 applied to the quadratic 
forms 

where m = n for generic ellipsoidal coordinates in E^c and m = n + 2 for the 
corresponding cyclidic coordinates. 

For coordinates corresponding to the direct sum of two spheres viz. (3.5) we can 
merely apply the result of the previous theorem for spherical coordinates to each 
of the pairs of quadratic forms 

p+i 2 p+i 

a=l 01=1 

n+2 

a=p+2 

Indeed the freedom to subject the coordinates uj and the Ci to birational transfor- 
mations in the expressions for generic cyclidic coordinates allows us to let ei — 00. 
(In this particular case the resulting coordinates can be identified with generic ellip- 
tical coordinates on the n sphere.) The process described in Theorem 1 enables one 
to pass from the elementary divisors to [Ni,N2, ....,Np],Ni > l,i = 1, ...,p, 

see [5,11]. It is then possible (via Theorem 2) to take ei — > 00 in which case 

[A^i, N2, Np] corresponds to the various generic coordinate systems in Euclidean 



n+2 2 



a=p+2 



U 
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consider the quadratic forms Cl, ^ corresponding to elementary divisors [21 • • • 1] 



VIZ. 

„2 



[A — ei)^ A — ex A — en+2 

Putting A — > 1 /A, Cj — > 1 /cj we find (with the use of ^ = 0) 



2 2 ii-r-^ I \ 



k 



or, in the hmit as ei — > oo, 

n+2 2 

A - Cfc 



fe=3 



Now if we perform the hmiting procedure (2.13) - (2.15) on L(u) (in n + 2 dimen- 
sions), let M — > 1/tt, Ci — > l/cj, let ei — > oo and then evaluate at pi = 0,a;i = 1 we 
find (up to a common factor u) 



n+2 

4 y 



Here, we have made use of the fact that J = 0. This agrees with (2.28) and with 
the L('u) operator given by Kuznetsov [4], and corresponds to elementary divisors 

oo 

[2,1,--- ,1]. 

This process can be generalised for quadratic forms corresponding to elementary 

oo 

divisors [iV, 1 • • • , 1], iV > 1. We have that 



n+2 9 N n+2 

xi 



h=u^ "^xi+u^ ^2xiX2+y^ xiXN-i- — - — , ^ = xiXN+i-i+ xi. 

i 

The corresponding L(u) operator is 



U — Ck 

k=N+l " i=l k=N+l 



n+2 

Sr 



where 



s\ = ^[n{0,p2, ...,Pn+2) - ^{i,X2, ...,Xn+2)], 

sl = ^[n{0,p2, ...,Pn+2) + ^{'i-,X2, ...,Xn+2)], 

1 ^ 
4= -2[^'i^^P2,-;Pn+2;'i-,X2,....,Xn+2)], S' = ^U^ ^ ^ XiPm- 

j=2 i+m=j 



16 E.G. KALNINS, V.B. KUZNETSOV AND WILLARD MILLER, JR. 

4. Branching Rules for the Construction of Orthogonal Non Generic 
Complex Integrable Systems on SnC and -E„c- To deal with the non generic 
separable coordinate systems in E^c and SnC we must combine coordinate systems 
for both manifolds [5]. (See [13,14] for tabulations of all cases for small values of 
n.) The branching laws for graphs on these manifolds are summarised below 

{Snc\ ••• 

(1) i 



PiC 



(2) i , A,>1 

Ep.c 

{E-ncl ••• l^il ■■■\ 

(3) i 



S. 



{End ■■■ |e: 



(4) i , A, >1. 

EpiC 

As an example consider the coordinate system given by the graph on S^ic, 

{S2C kil 65] 

i 

[E2C 1/3 1 /4I 

The coordinates for this graph can be obtained from those of the generic graph 

('S'4C lexical 64 1 65 1 

via the limiting transformations 



Cj = ei + e + e^fj, j = 3, 4, /s 7^ / 



4 



Uj = ei + e + e Uj, j = 3,4: 



where e 0. 

The corresponding coordinates are then given implicitly by 

_ 2 _ (^1 - ei){u2 - ei) 2 _ (^1 - 65) (^2 - 65) 

65-61 (ei-65)^ 

(wi-ei) {u2-ei) (-ui - 6i)(w2 - ei) 

-2xiX2 = 1 7 

65-61 65-61 (65-61)"= 

{ui - ei){u2 - ei) _ _ 

-[f3 + f4-U3-U4) 



x^ 



65 - ei 

2 {ui-ei){u2-ei) {u3- fs){u4- fs) 



(es - ei) if 3 - U) 



_ 2 _ (m - ei)(tt2 - ei) {U3 - f4){u4 - fi) 

Xa / s 



/I* I* \ 
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Using the nomenclature given previously we have that L(w) has the two different 
forms 

Li(w) = 1 1 , U = Ul,U2 

[u — eiY u — ei u — ei u — 
L2(it)=ZiH — H — , u = us,U4. 

As usual, the separation variables are the zeros of the equation L\{u) +iL\{u) = 
with u = ui,U2 when A = 1 and u = us,U4 when A = 2. Each of the L;^ separately 
satisfy (2.12) for A = 1, 2. In addition we have 

{Li.Li(«),L2-L2(^)} = 0. 

This example illustrates how to derive the substitutions that enable the various 
branching laws to be obtained from a generic form. For the sphere SnC and generic 
coordinates 

(S'nc|---|eJ^^~*|ei| • • • |es| • • • I q> s, X>1, 
if we make the substitutions 

efe = eo + e«-^+' + e^+'^'^fk, k=l---s, 

Uj = eo + e«-*+^ + e^+'i-'uj, j = l---q 

where e — > 0, then the graph illustrated transforms into 

{Snc\ |eol 

i 

[^.Cl fl\ \fs\ 

The remaining branching rule is obtained from a graph of the type 

(5'nc|---|eo|ei|....|es|. 

By means of the substitutions 

Cj = eo +e(/o - /j), j = Uk = cq + e{uk - fo) , A; = !,••• ,p 

we obtain the coordinate system coming from the graph 

{Snc\ |eo| 

i 

{Spc\ fo\ 

The corresponding substitutions for the analogous Euclidean space branching 
rules are essentially identical. To completely specify the coordinate systems on 
these manifolds we need a few more substitution rules. Firstly consider the graph 

[Epc\hi\....\hp\ + [£;,c|/i|....|/,|. 
This graph can be obtained from the generic graph for £^(p+q)c via the substitutions 

Ui = \-Ui, Cj = \- hi, i = l,-- ■ ,p 

e e 

K2 _ K2 
Ui = — + Ui, ei = — + fi-p, i = p+l,--- ,p + q, Ki^ K2. 

The graph [£^qc|^i|----|^g|) q < n can be obtained from [£^„c|^i|----|^n| via the 
substitutions 

Ui = ei + e'^~'^~^ + e'^~'^Ui, z = 1, • • • , n 

= ei + e^"''"^ + e^-'^/ij, j = 1, • • • , 5 ek = hk, /c = g + 1, • ■ ■ , n. 

Given these substitution rules all the corresponding graphs for E^c and can 
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5. Quantum Integrable Systems on Complex Constant Curvature Spaces 
and the Quantum Gaudin Magnet. To deal with the quantum version of this 
description of separation of variables we consider the Schrodinger or Helmholtz 
equation [3]. 

where E for the moment. Separation of variables means (roughly) the solution 
of this equation of the form 

(5.2) * = n^=iV;«(ya;/ii,--- ,M 

where the quantum numbers hj are the eigenvalues of mutually commuting opera- 
tors 

[A, Aj] =0, z, j = 1, • • • , n, A = 7^, Aj^ = hj^. 

Furthermore these operators can be represented as symmetric quadratic elements 
in the enveloping algebra of the symmetry algebra of the equation (5.1) in the case 
of SnC and Enc- The standardised representations of these symmetry algebras are 

SO{n + 1) : Map = XaPp - XpPa, Pa = a, P = 1, ■ ■ ■ , u + 1. 

(5.4) [Map, Mjs] = Sa^Msp + SasMp^ + Sp^MaS + SpsMja 

d 

E{n) : Map, ^i = -Q^ 

(5.5) [Map, V^] = Sp^Va - SajVp, [Va,Vp] = 

where [, ] is the commutator bracket. Much of the analysis goes through as it did 
in the classical case with, of course, some critical differences. For the quantum 
Gaudin magnet one considers the sum of rank 1 Lie algebras A = ©^=iS0q,(3) 
where the generators of the algebra satisfy the commutation relations (1.6) and the 
inner product is defined as in §1. The Casimir elements of A have the form 

(5.6) {Sa,Sa) = ka{ka + l) 

where ka is a constant when the generators of sOq:(3) determine an irreducible rep- 
resentation, [12]. The quantum Gaudin magnet is the quantum integrable Hamil- 
tonian system on A given by m commuting integrals of motion : 

Ha =2 ^ [Ha, Hp] = 0. 

This is the m site SO{3,C)-XXX quantum Gaudin magnet, [3]. These integrals 
(operators) satisfy 

mm m 

(5.7) J2'^a = Q, Xl^«^« = + 



SEPARATION OF VARIABLES AND THE GAUDIN MAGNET 



19 



where J = Sq, is the total momentum operator. The complete set of com- 

muting operators consists of Ha,^^ and J^. The integrals are generated by the 
2x2 operator L{u) given by (1.11) understood in the operator sense. The quantum 
determinant is 

(5.8) q-detL(«) = -Aiuf - '-{Biu), Ciu)} = - - ^f^^^^ 

a=l a=l ' 

with L-operator satisfying the r matrix algebra 

(5.9) {L{u)X{v)\ = -^{pX{n)^i{v)l i = 

u — V 

(compare with (1.13)). The total momentum J has components 

(5.10) JT^ = ^(p-p + x-x), JT' = ^(p-p-x-x), J-^ = ^(p-x + x-p) 

The Lie symmetries and J form the direct sum so{m) ® so(3) with com- 

mutation relations 

(5.11) [Map, X\ = 0, \J\ J^] = iejuJ^ 
subject to the constraints 

(5.12) 3^ = lY,Mlp + ^m{m-A), 2( + z J^) = x • x. 

a</3 

Considering again separation of variables using the coordinates of the irreducible 

block (2.1), the separation variables are defined, as before, as the zeros of the off 
diagonal elements B{u) of the L matrix. The g-determinant is the generating 
function of the commuting integrals of motion 

(5.13) [q-detL(w), q-detL(v)] = 0. 

On the n sphere the algebra A is realised by taking the canonical operators 

(5.14) 4 = ^(^a +^a): 4 = ^(Pa-^a): 4 = ^{Pa, ^a}, 1, 1 

where (Sq,,Sq,) = Furthermore 

(5.15) {^o.,sp)=^-{Mlp + \) 

and this establishes the relationship between the Sq, of A and the Aiap of so(n + 1). 
The integrals transform into the family of integrals 

(5.16) ^ = ^ KUa = \ E -^^{Ml, + h 
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n+1 2 
X7 



and the coordinates are given by 

For each uj variable there is defined the conjugate variable (operator) 

1 '^"^^ 1 

(5.17) Vj = -iA{Uj) = jYl {Pa,Xa}- 



a=l 

Separation of variables is then the process of changing to the new variables f j, Wj, c 
and . In fact we have that 

(8.18) ltf„^„HfHJ3_g_l_^^,„^j, 

The separation equations can be obtained by substituting u = Uj into the q-detL(u), 
making the choice 

and looking for the solutions of the spectral problem Tia^ — . The separation 
equations are then 

where 

* = n^+V,(«,), h'^-h^ + lj^T^- 

These are the separation equations for the Helmholtz equation in generic ellip- 
soidal coordinates. Note that for this choice of Vj we have taken (sq,, Sq-) = —3/16. 
The L operator is given in direct analogy with the classical case, 

n+l 



(5.21) Uu) = Y.77Z 



a.=\ 



The treatment of generic ellipsoidal coordinates on E^c follows similar lines, 
with the L operator given in this case by 

The limiting procedures described in the classical case work also in the quantum 
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In the case of cyclidic coordinates we can adopt the same strategy as in §3: we 
impose the conditions J" = and proceed as before. The natural setting in this case 
is again to use hyperspherical coordinates. The total momentum has components as 
in (5.10). If Pol, Xa are now vectors in hyperspherical coordinates then we can derive 
the standard quantum r matrix algebra as above with the same formulas valid. The 
constraints are now 2{J'^ + ij"^) = x • x = 0, = |^(^ — 4). Coordinates Uj and 
their conjugate operators Vj can be chosen as before. In particular if we make the 
choice 

Q ^ n+2 ^ 

, , _ a 1 1 

J a=l ■' 

and look for solutions of = of the form * = (j^'^~'^^l'^Wy^y'^i{ui), then we 
obtain the separation equations which coincide with the the equations given by 
Bocher,[ll], viz. 

n+2 



3 a=l J ^ J 



^{n' - 4)uJ - ^(2n - n'){^ltl e.)u]'' + E^o ^P^^ 



^3 



for suitable 7^. Note that the solutions ^ arc not strictly separable in this case 
but are what is termed R separable, i.e., separable to within a non separable factor 

^(n-2)/2^ [1]. 

We conclude this work by pointing out that one can obtain a complete set of 
constants of the motion associated with an orthogonal separable coordinate system 
{wi, ■ ■ - Un} for the Schrodinger equation (5.1) directly from (5.1) itself, [15, 16]. It 
is well known that all orthogonal separable coordinate systems for the Schrodinger 
equation on a Riemannian space are obtainable via the Stackel construction, e.g., 
[17, 18]. Thus if {u} is a separable orthogonal coordinate system for (5.1) there 
exists an 71 X n nonsingular matrix S = {S~^_p{ua)) such that du S^p = if 7 7^ ct, 
and such that the nonzero components of the contravariant metric tensor in the 
coordinates {u} are gf""(u) = T"^"(u), a = 1, • • • , n, where T is the inverse matrix 
to S, 

n 

(5.24) 5]T"^^(u)^^^(«^) = 5«. 
The constants of the motion are 

n 

(5.25) ^^ = ^T/3«(u)(a2^+/«(««)5„J, An = 'H. 

01=1 

Here, /«(««) = ln(V^5r««). (The fact that 
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follows for any space of constant curvature by noting that it is equivalent to 
the statement that the off diagonal elements of the Ricci tensor must vanish for 
an orthogonal coordinate system: R^fs = 0, a ^ f3, [17].) One can show that 
[^Q,^/?] = 0, [17, 18]. Furthermore, if = Il'^^iipa{ya',hi,- ■ ■ ,hn) satisfies the 
separation equations 



n 



(5.26) a2^^' + /«(w«)a„„^' = 5^5«/3(««)/i^*, a-1,- 
then 

(5.27) Ai3'^ = hfi^, /3=l,---,n. 

Now fix 7, 1 < 7 < n, and denote by y the coordinate choice 

y = (ui, • • • ,u^_i,r, -u^+i, ■ ■ ■ ,'Un) 

where r is a parameter. We see from (5.26) that if ^ satisfies the separation 
equations then 

n 

(5.28) V^;' {dl_^^^ + Uu^)du^,l;^) U,=r^'(u) = S^p{T)Ap^{u). 

13=1 

On the other hand, from (5.1) we have 

^^-7' {dl^i^i + Uu,)du,^^) U,=r*(u) = -J^(7i(u)-7i^(y))*(u) 
where 

This suggests the operator identity 

1 " 

(5.29) -— - {n{n) - H\y)) = S,p{T)Ap, 

9 yy) 

i.e., that the left-hand side of (5.29) is a one-parameter family of constants of the 
motion and that as r runs over a range of values and 7 = 1, • • • , n the full space of 
constants of the motion associated with this separable system is spanned, [15, 16]. 
Now (5.29) is equivalent to the conditions 

and these conditions are easily seen to follow from (5.24) and from 

n 

Similarly, for the classical case an expression for the Hamiltonian analogous to 
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